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Chiral dynamics of the nuclear equation of state 
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We present a new chiral power expansion scheme for the nuclear equation of state. The 
scheme is effective in the sense that it is constructed to work around nuclear saturation 
density. The leading and subleading terms are evaluated and are shown to provide an 
excellent equation of state. As a further application we considered the chiral quark con- 
densate in nuclear matter. Already at nuclear saturation density we predict a substantially 
smaller reduction of the condensate as compared to conventional approaches. 
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1. Introduction 

00 . 

The nuclear equation of state lies at the heart of nuclear physics. It is therefore of 
considerable importance to find the appropriate form of chiral perturbation theory (%PT), 
the most powerful tool of modern nuclear physics, and apply it systematically to the 
' nuclear many body problem. 

A key element of any microscopic theory for the nuclear equation of state is the el- 
>• . ementary nucleon nucleon scattering process. In the context of this problem was 
first addressed by Weinberg who proposed to derive a chiral nucleon nucleon potential in 
time ordered perturbation theory Jl|. The NN-phase shifts then follow from the solution 
of the Schroedinger equation appropriately fed with the chiral potential. This program 
was carried out by Ordonez, Ray and van Kolck 0. Weinberg's scheme is plagued by 
two problems. First the systematic renormalization of the chiral potential scheme is an 
open problem and not resolved. The scheme is regularization scheme dependent. Sec- 
ond the use of a potential approach in nucleon nucleon scattering can be justified only if 
retardation effects are small. 

To overcome both problems the present author proposed first to apply chiral power 
counting rules to the nucleon nucleon scattering amplitude directly [[|. This approach 
cures both problems. First the manifest covariant version of the chiral Lagrangian can 
be applied. A relativistic form of was applied before by Gasser [[| in the 1-nucleon 
sector but rejected and replaced by the heavy mass formulation of ||. There are 
two problems inherent with the relativistic approach of [ffl. First, any covariant derivative 
acting on the nucleon field produces the large nucleon mass, m^, and therefore must be 
assigned the minimal chiral power zero. Thus an infinite tower of interaction terms need 
to be evaluated at a given finite chiral order. Second, the straightforward evaluation 
of relativistic diagrams involving nucleon propagators generates positive powers of the 
nucleon mass from loop momenta larger than the nucleon mass. The chiral power counting 
rules are spoiled. 
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The two problems are solved as follows. Rather than performing the 1/rriN expansion 
at the level of the Lagrangian density, as achieved in the heavy mass formulation of 
%PT, it is equally well possible to work out this expansion explicitly at the level of 
any individual relativistic Feynman diagrams [[J. The second problem is overcome by 
a suitable reorganization of the infinite tower of chiral interaction terms. Of course we 
expect our relativistic scheme to be equivalent to the heavy fermion formulation of chiral 
perturbation theory in the one-nucleon sector. 

The second important observation made first in || is the fact that the chiral power 
counting rules can be generalized for 2-nucleon reducible diagrams. For a given diagram 
each pair of intermediate nucleons causes a reduction of one chiral power as compared to 
the 'naive' chiral power. The non perturbative structures like the deuteron bound state 
are generated naturally since the bare 2-nucleon vertex is renormalized strongly so that 
it effectively carries chiral power minus one. 

An attempt to apply the generalized chiral power counting rules of to the nuclear 
many body problem quickly reveals that even though the pion dynamics remains pertur- 
bative the local 2-nucleon interaction requires extensive resummations. 



2. Chiral expansion scheme for nuclear matter 

It is straightforward to combine the chiral expansion of the nucleon scattering ampli- 
tudes with the low density expansion. Naive application to the energy per particle, E{kp), 
of isospin symmetric nuclear matter would lead to an expansion of the form: 

«M - ? *® (!)" U> 

with A ^ m n and p = 2k F /(3ir 2 ). The expansion coefficients E n are functions of the 
ratio kp/m n . Such a scheme is obviously restricted to extremely low density, if useful at 
all, since the typical scale A is small. This follows for example if one considers the 2- 
nucleon rescattering contribution to the 3-nucleon scattering process. Then the deuteron 
pole term induces structures like k F j \fm^ Id where €p> — 2 MeV denotes the deuteron 
binding energy. One is lead to identify A ~ ^m N e D . 

A scheme working at nuclear saturation density therefore requires a further resumma- 
tion. Obviously an expansion of the form 

**> - ®" <» 

with typical small scales, As, like y/m,N en and typical large scales A^ ~ 4:71 f n ~ 1 GeV 
must be achieved. The expansion coefficients E n are complicated and hitherto unknown 
functions of the Fermi momentum kp. They can be computed in terms of the free space 
chiral Lagrangian furnished with a systematic resummation technique |5]]. 

In this work we present a somewhat less microscopic approach in spirit close to the 
Brueckner scheme but more systematic in the sense of effective field theory. Since the 
typical small scale As is much smaller than the Fermi momentum kp ~ 265 MeV at 
nuclear saturation density, one may expand the coefficient functions E n around kp in the 
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following manner 
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Note that we do not expand in the ratio m^jkp. If one expanded also in this ratio m n /k F 
one would arrive at the Skyrme phenomenology ]7| applied successfully to nuclear physics 
many years ago. It should be clear that this scheme is constructed to work around nuclear 
saturation density but will fail at small density. We note that also conventional approaches 
like the Walecka mean field or the Brueckner scheme are known to be incorrect at small 
density. 

Technically our scheme can be generated by the effective Lagrangian density 
Cint(k F ) = |f iV 75 7 m -(^tt) -tN 

+ (9o(k F ) + \g 2 A )(Nl 5 T 2 C- 1 N t )(N t CT 2l5 N) 

+ ^ (g 1 (k F ) + ~g^(N^TT 2 C- 1 N t )(N t CT 2 r^N) (4) 

where the couplings go = go(k F ) and g\ = gi(k F ) are density dependent. A more sys- 
tematic derivation of the expansion (||]) and (ff) applying suitable resummation techniques 
will be presented elsewhere 0. 

The chiral power counting rules are simplified significantly as compared to a fully mi- 
croscopic scheme. The presence of a further small scale k F ~ Q ~ m n does not anymore 
generate an infinite tower of diagrams, to be considered at a given chiral order, since 
by construction the troublesome local 2-nucleon vertex must not be iterated, i.e. terms 
proportional to g^ikp) and g^ikp) with n > 1 are already included in go(kp) and g±(k F ) 
and therefore must not be considered. The pion dynamics, if properly renormalized, re- 
mains perturbative like in the vacuum case. Consider for example the two loop diagrams 
depicted in Fig. |l| where the nucleon line with a 'cross' represents the projector onto the 
Fermi sphere 

ASW(p) = (l-p + m N )27TtQ(p )5(p 2 -m 2 N )e(k 2 F -p 2 ) . (5) 

The first diagram in Fig. [I] is proportional to go y i(k F ) k F and is therefore ascribed the 
chiral order Q 6 since the effective vertex (?o,i(^f) ~ Q° carries chiral power zero. The 




Figure 1. Leading contribution of chiral order Q 6 . 
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Figure 2. Subleading contribution of chiral order Q . 

second diagram, the one pion exchange contribution, is also of chiral order Q 6 since it 
is proportional to k F multiplied with some dimensionless function /(m^/fci?)- In Fig. ^| 
we collected all diagrams of chiral order Q 7 . Here we introduced two types of 2-nucleon 
vertices. The filled circle represents the full vertex of (|J) proportional to g> ,i + 9a/^ 
and the open circle the counter term proportional to g\/&- The dashed line is the pion 
propagator and the directed solid line the free space nucleon propagator. We point out 
that the diagrams b), d), f) and h) in Fig. |2| are divergent. The leading chiral contribution 
of the sum of all diagrams, however, is finite. 

It is instructive to also display the non relativistic form of the Lagrangian density (^) 
as suggested by Weinberg |j] in the context of the nucleon nucleon scattering problem 

dnt{k F ) = ^N(a-V)U-f)N 

+ ^Jt (go(k F ) + - g\j (N a a 2 r 2 N l ) (iV* r 2 a 2 a N) 

+ ^ (gi(kF) + \g 2 A ) (Na 2 tt 2 N 1 ) (N'^r^N) (6) 

with now a two component spinor field and N = iV" for notational convenience. Ac- 
cording to Weinberg the two particle reducible diagrams are to be evaluated with static 
pions. We therefore evaluate all diagrams of Fig. |^ with the interaction vertices of @ and 
static pion propagators. The solid line with a 'cross' now represents the non relativistic 
limit of (P) and the directed solid line the free non relativistic nucleon propagator. We 
remind the reader that in the relativistic approach (^) we perform the chiral expansion 
at hand of a given Feynman diagram. For technical details we refer to || . Both schemes 
indeed lead to identical results for all diagrams except diagram h). Only in the extreme 
low density limit we find agreement of the non relativistic with the relativistic scheme 
also for diagram h). We conclude that it is incorrect to use static pions in two-particle 



5 



t ' r 



static ji's 




i i i i i i i 

0.0 0.5 1.0 1.5 2.0 



P [P ] 



Figure 3. The equation of state for isospin symmetric nuclear matter. 

reducible diagrams. Therefore the use of a potential in the xPT approach of the nuclear 
many body problem or the nucleon nucleon scattering problem cannot be justified. A 
relativistic approach is required. 

In Fig. |] our result for the isospin symmetric nuclear equation of state is shown. We 
compare our relativistic approach (f|) with the static potential approach of @. The 
relevant coupling g + g 1 is adjusted to obtain nuclear saturation at k F = 265 MeV. The 
non relativistic scheme gives a rather poor equation of state as shown in Fig. [|. The 
binding energy and the incompressibility are too small. We emphasize that the coupling 
functions <7o,i(&f) are to be determined from the nuclear equation of state. It is therefore 
legitimate to cure the 'non relativistic' equation of state by giving the coupling a residual 
density dependence 

*<*,)+»<*,) = 7o + o4(l-Q + 0,(|) 2 (l-Q 2 + .... (7, 

However there is a strong consistency constraint: according to our scale argument (|3]) 
the density dependence of the couplings <?o,i(&f) must be weak. If the nuclear saturation 
required a strong density dependence our scheme had to be rejected. We turn to the 
relativistic scheme. Here the set of parameters go+gi — 2.7, g& ~ 1.2, m n ~ 140 MeV and 
f n ~ 93 MeV give an excellent result for the equation of state. The empirical saturation 
density kp ~ 265 MeV and the empirical binding energy of 16 MeV are reproduced. The 
incompressibility with k ~ 234 MeV is also compatible with the empirical value (210 ±30) 
MeV of ||. We conclude that in fact the consistency constraint points towards the correct 
scheme @. 
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3. The chiral order parameter (qq) 

The quark condensate, (qq)(p), is an object of utmost interest. It measures the degree 
of chiral symmetry restoration in nuclear matter. Furthermore it is an important input for 



QCD sum rules or the Brown Rho scaling hypothesis [10]. According to the Feynman 
Hellman theorem the quark condensate can be extracted unambiguously from the total 
energy, E(p), of nuclear matter once the current quark mass dependence of E(p,i71q) is 
known 

(<?<?> (P) = ^T^-E(p,m Q ) (8) 

with niQ = m u = nid- Recall that E(p)/V = (mjv + E(p)) p is determined by the nuclear 
equation of state E(p). Since our chiral approach was set up to treat the pion dynamics 
and therewith the current quark mass dependence of the equation of state properly it is 
very much tailored to be applied to the quark condensate. 

It is convenient to consider the relative change of the quark condensate since it is 
renormalization group invariant: 

(qq)(p) = j _ ^Np _ aMp (q) 
(qq)(0) mlfl 2m* ft ' U 

The second term in (|9|) follows from the nucleon rest mass contribution to the total energy 
of nuclear matter together with the definition of the pion nucleon sigma term, X^: 

Z«N = m Q (N\qq\N) = m Q ^ . (10) 

drriQ 

The last term in (Q) measures the sensitivity of the nuclear equation of state, E(p,m n ), 
on the pion mass 

<*M = -^^j^E(p,m w ) = (l + 0(ml))-^-E(p,rn 7T ) (11) 
{qq){0) orriQ v v // Qm^ 

where we consider now the current quark mass tuq = mQ(m n ) as a function of the physical 
pion mass. We emphasize that the second term in (|9]) written down first in [|IT1 -|T3|| does 
not probe the nuclear many body problem and therefore should be considered with great 
caution. It is far from obvious that this term is the most important one at nuclear 
saturation density. The only term susceptible to nuclear dynamics is the last term in (|S|). 

In Fig. H we present our result for the quark condensate in nuclear matter. We confront 
the 'leading' term driven by E^-tv — 45 MeV with first, our result applying static pions 
and second, our result including important retardation effects. The former result is rather 
close to the 'leading' order result. The nuclear many body dynamics as induced by static 
pions appear to effect the quark condensate little. This confirms results obtained in 
the Brueckner [M| and Dirac-Brueckner [15] approach. Note that neither the Brueckner 



nor the the Dirac-Brueckner approach includes pionic retardation effects systematically. 
We celebrate the central result of this section: as shown in Fig. ^ the inclusion of pionic 
retardation effects leads to a significantly less reduced quark condensate in nuclear matter. 

Our result implies strong consequences for the QCD sum rule approach to the properties 
of hadrons in nuclear matter. 
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Figure 4. The quark condensate in isospin symmetric nuclear matter. 
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